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Abstract 

Data from the literature were used to calculate 
separate regressions of summer chlorophyll a (Chla) 
concentration on spring total phosphorus (TP) concentration 
for deep lakes (i.e., lakes which remain thermally 
stratified during the summer) and shallow lakes (i.e., lakes 
which mix intermittently during the summer). Significant 
differences were found in the spring [TP]-summer [Chla] 
relationship for the two lake types (P < 0.05). The mean 
ratio of summer [TP] to spring [TP] was also significantly 
different in deep and shallow lakes (P < 0.001). This 
difference in the summer [TP] to spring [TP] ratio is the 
explanation offered for why the spring [TP]-summer [Chla] 
relationship was different in the two lake types. Internal 
loading of phosphorus from the sediments to the surface 
water in Shallow lakes was suggested as the cause of the 
difference in the summer [TP] to spring [TP] ratio in the 
two lake types. 

The contribution of internal loading from the sediments 
to the TP budget and to the [TP] in the surface water was 
investigated from May to November in Nakamun and Halfmoon 
Lakes, two shallow, productive lakes in Alberta. During the 
Summer, Nakamun Lake was intermittently stratified and 
Halfmoon Lake was weakly, thermally stratified. While the 
lakes were stratified, the water overlying the sediments was 
anoxic and TP levels increased in the deep water. During the 


stratified periods, the sediments were the main source of 
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TP, contributing 1468 and 147 kg of TP to Nakamun and 
Halfmoon Lakes, respectively, During these same periods, 
external loading to Nakamun and Halfmoon Lakes was 
insignificant (37 and 5 kg, respectively). Average release 
rates of TP calculated from the in-lake TP budgets for 
Nakamun and Halfmoon Lakes (12.7, 15.6 mg/m*/day, 
respectively) were slightly higher than the rates predicted 
from laboratory experiments on sediment cores (9.7 and 7.2 
mg/m?/day, respectively). Quantitative estimates of vertical 
water exchange in Nakamun Lake and in-lake TP budgets for 
both lakes indicated that phosphorus released from the 
sediments was transported to the surface water during 
mixing. After eight of the nine mixing events which 
immediately followed stratified periods, the [TP] increased 
3-43% and 31-52% in the surface water of Nakamun and 


Halfmoon Lakes, respectively. 
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Figure 2. Mean total phosphorus concentration (+ SE) of the 
water over lake sediment cores versus time. Cores were 
incubated in the dark under: (A) anoxic and oxic 
conditions at 10° C and (B) under anoxic then oxic 
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Figure 3. Mass of total phosphorus above and below the depth 
of 2.5 m versus time in: (A) Nakamun and (B) Halfmoon 
Lakes. Mass of total dissolved phosphorus (TDP) below 
2.5 m also given for Nakamun Lake. Black bars indicate 
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Figure 4. Three dimensional graphs of time versus depth 
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Figure 5. Three dimensional graphs of time versus depth 
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Figure 6. Total phosphorus concentration ([TP]) and 
chlorophyll a concentration ([Chla]) in the trophogenic 
zone from May until October, 1982 in: (A) Nakamun 
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Figure 7. Percent composition of dominant groups in the 
algal community in the trophogenic zone of Nakamun Lake 
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I, Introduction 
The nutrient phosphorus (P) is the primary factor which 
limits summer algal biomass in north temperate lakes 
(Schindler 1977). Consequently, total phosphorus (TP), often 
measured at spring overturn, iS a good predictor of summer 
chlorophyll a (Chla) concentrations (Sakamoto 1966; Dillon 
and Rigler 1974; Oglesby and Schaffner 1978; Smith 1982; 
Prepas and Trew 1983). There are suggestions that the summer 
[Chla] to spring [TP] ratio is higher in shallow lakes 
(i.e., lakes which mix intermittently during the summer) 
than in deep lakes (i.e., lakes which remain thermally 
Stratified during the summer) (Oglesby and Schaffner 1975; 
Schaffner and Oglesby 1978; Prepas and Trew 1983). 

The proposed explanation for the higher summer [Chla] 
per unit of spring [TP] in shallow lakes is that TP 
concentrations increase from spring to summer in the surface 
water of shallow, but not deep, lakes (Prepas and Trew 
1983). External loading cannot account for the observed TP 
increases in many shallow lakes during summer (Larsen et 
al. 1975; Stevens and Gibson 1977; Reynoldson and Hamilton 
1982). Internal P loading is often considered responsible 
for these increases, uSually through P release from 
profundal sediments when the overlying water 1S anoxic, with 
subsequent transport of this P to the surface water (Larsen 
et al. 1981: Stefan and Hanson 1981). 

Alberta has many shallow, productive lakes. In these 


lakes, TP concentrations in the surface waters increase 
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throughout the summer (Prepas 1983; Prepas and Trew 1983). 
Prelimnary studies on some of these lakes suggest release of 
phosphorus from the sediments under anoxic conditions is the 
most significant source of internal phosphorus loading 
(Prepas and Wisheu unpublished data). 

In this study, I used data gathered from various 
sources to compare the spring [TP]-summer [Chla] and summer 
[TP]-spring [TP] relationships in deep and shallow lakes. It 
has been suggested these relationships are different in the 
two lake types, but these differences have not been tested 
Statistically. I also examined the P loading on two shallow, 
productive lakes in Alberta and conducted laboratory studies 
on P release from the sediments to determine the importance 
of internal loading from the sediments on the TP budget and 
[TP] in the surface water of shallow lakes. This is the 
first attempt to quantify the input of P release from the 


sediments in Albertan lakes. 
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II. A Comparison of the Phosphorus-Chlorophyll Relationships 


in Shallow and Deep Lakes 


A. Introduction 

The relationship between the total phosphorus (TP) and 
chlorophyll a (Chla) concentrations in lakes has been 
examined by a number of authors (Sakamoto 1966; Dillon and 
Rigler 1974; Oglesby and Schaffner 1978; Smith 1982; Prepas 
and Trew 1983). This relationship is fairly consistent when 
spring or winter [TP] is used to predict the average summer 
[Chla] for a lake (Nicholls and Dillon 1978). The spring 
[TP]-summer [Chla] relationship has been used as a simple, 
yet valuable tool in lake management (Dillon and Rigler 
(97 50% 

The applicability of models to predict summer [Chla] 
(or some other estimate of the phytoplankton standing crop) 
would be improved if the residual variance in these models 
were reduced (Nicholls and Dillon 1978: Smith 1982). There 
are suggestions that the variation in the spring [TP]-summer 
[Chla] model would be reduced if the lakes were divided into 
two categories: deep lakes (i.e., lakes which remained 
thermally stratified during the summer) and shallow lakes 
(i.e., lakes which mix intermittently during the summer) 
(Oglesby and Schaffner 1975; Schaffner and Oglesby 1978; 
Prepas and Trew 1983). 

Further analysis of Prepas and Trew's data shows that 


the summer [Chla] to spring [TP] ratio is higher in shallow 
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lakes than in deep lakes. The measured Chla concentrations 
in six of the eight shallow lakes in their study were 
greater than the values predicted by Dillon and Rigler's 
(1974) spring [TP]-summer [Chla] model. The measured summer 
Chla concentrations in only 7 of the 18 deep lakes in their 
study were greater than the predicted values. 

Oglesby and Schaffner (1975) calculated separate linear 
regressions of summer [Chla] on winter [TP] for deep and 
Shallow lakes. They assumed that the spring [TP]-summer 
[Chla] relationships were different in the two lake types. 
Oglesby and Schaffner's regression line for shallow lakes 
has a greater slope and intercept than the regression line 
for the deep lakes. However, their data sets were too small 
to test whether these differences were significant. 

During summer, the [TP] in the trophogenic zone 
increased in shallow lakes and decreased in deep lakes in 
Prepas and Trew's study. Differences in summer [TP] to 
spring [TP] ratios in the two lake types could explain why 
Shallow lakes often have higher summer [Chla] than deep 
lakes, relative to spring [TP]. 

To examine whether the differences in the spring 
[TP]-summer [Chla] relationships in deep and shallow lakes 
were Significant, I calculated separate regression lines for 
summer [Chla] on spring [TP] for both lake types. I then 
compared the regression lines for deep and shallow lakes 
with each other and with Dillon and Rigler's (1974) 


regression line by analysis of covariance. I treated Dillon 
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and Rigler's regression line as a model for deep lakes 
because most of the lakes in their data set were deep. The 
results support the hypothesis that the spring [TP]-summer 
[Chla] relationships are different in deep and shallow 
lakes. Differences in the summer [TP] to spring [TP] ratios 
in the two lake types offer the best explanation for the 


differences in the spring [TP]-summer [Chla] relationship. 


B. Methods and Materials 
The data for the regression lines were gathered from 
several sources: Dillon and Rigler (1974), Oglesby (1977), 
Oglesby and Schaffner (1978), Prepas and Trew (1983), C.E. 
Gibson (University of New Ulster, Coleraine, Northern 
Ireland, pers. comm.), and E. Mills (Cornell Biological 
Field Station, Bridgeport, N.Y., pers. comm.) (Table 1). The 
data consisted of: 
ky) SpeingiftTPliives! theelake [TP] eabespringioverturnt 
Phosphorus data collected in late winter (Oglesby and 
Schaffner 1978) were treated as spring values since they are 
comparable (Nicholls and Dillon 1978). 
2) [Chla] i.e., average lake concentration of chlorophyll a 
during the summer. 
3) Summer [TP] i.e., the average summer [TP] in the 
trophogenic zone. The trophogenic zone was defined as the 
water above the depth of 1% surface-penetrating irradiance. 
I did not use some data from the original data sets in 


our analyses. Mills (pers. comm.) and Oglesby and Schaffner 
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(1978) both provided data for Oneida Lake in 1975. Mill's 
data were used for the calculations because they were based 
on a much larger sample size than Oglesby and Schaffner's 
(23 versus 2 samples for spring [TP] and 65 versus 2 samples 
for [Chla], respectively). Three lakes from Prepas and 
Trew's data set were excluded from the calculations. The 
spring [TP] values in two lake's (Joseph and Cooking) were 
well beyond the range of the other lakes in Table 1 (706 and 
303 mg/m*, respectively). A third lake, Hasse, was excluded 
because it could not be clearly classified as either a deep 
or shallow lake. | 

The data were transformed to log;, before regression 
analyses as in other spring [TP]-summer [Chla] models 
(Dillon and Rigler 1974; Prepas and Trew 1973). The 
regression analyses were performed with BMDP statistical 
packages (Dixon 1981) on the Amdahl model 580/5860 computer 
at the University of Alberta. The other statistical analyses 
were from Zar (1974). 

The data in Table 1 were used to calculate the 
regressions of summer [Chla] on spring [TP] for the two lake 
types. A regression line based on Dillon and Rigler's (1974) 
data set was also included in the analysis of covariance. 
The spring TP concentrations in 5 shallow lakes were higher 
than for any deep lakes. The spring [TP] in eight deep lakes 
were below the range of values for shallow lakes (Table 1). 
In analysis of covariance the range of the independent 


variaables shoulaybehsimigaretor aliidata|)sets; Thus lakes 


beasd 978M 
e'rant ede’ bam gdestoo. nu 
aslamsa’ $ steele Ter 
Bite enensit ters eda a Ra . s 
od? .nnoltaluqhen 26% so re occas we seas 7 
crew (pattoed Brantgaman) 9 'edet Oe) ma eoulav (80) ¢ * oe 
has 20°) | sidwh nd giedel sedge ef? Yo eonss edz Bnoyed Iie 2 
bebuioxe esw sae ~otel Boidt «A .(ylevideogwes «ta\om aa <4 
qaeb o ssdgis 26 be? tiapel> ¢iapals: meates Bivao of seuensd eu | 
Ps oe asi volieda 7 
 golegetpes sibted ‘oreo od Beaiatenssy. elev sieb.edtT > nm | 
eishon (91d0] dvemowe-f4¥] god 2q2: 1]djor nies S0eTIAR 
ad? (ENS wor? Sap Baqest iter aelgin bre MOLLE). 
[eoisacvese SOMA). aaiw bebroloeg: 2x8 sse¢isan' notadexped 
tesugmes 0962088 lebom ivetmA sdgino (° 8Or aoxid) caged ag on 
eseyisne Lesiveldsga sedge GAT sagtedlA io ySisievial oft 6. Fs 
‘ Los i /) ,(E°TEr) tek moTt erst OP 


ed3 sielyoties o7 enh <aias | sidst nk ‘e3eb efT 4 

ois! ows off FOP ESI onttae no) bela) stot io snpleasiger: 

(ave) a rehols bos aol lid ao Beead onit noiaerges! A “saat. 
-ooaizeves to eteylsas edd ai Sabutoni cela. ean sons : 

sedpid 9x69 een oe * pad . 

' aoa OSE 8 

(1 eisie?) ax 


; a vaphers ; 
neat ea a 


Re yy 


ae 
7 oi ; 
na ai 


- 


with more than than 10 mg/m? or less than 100 mg/m? of 
spring [TP] were excluded from the analysis of covariance. A 
range of 10 to 100 mg/m? of spring [TP] was chosen because 
it did not reduce the size of the data sets too much and 
there were no extreme outliers within this range. 

To compare the relationship between summer [TP] and 
spring [TP] in the trophogenic zone, and summer total 
nitrogen concentration ([{TN]) and summer [TP] in the 
trophogenic zone of the two lake types, I calculated the 
summer [TP] to spring [TP] and summer [TN] to summer [TP] 
ratios for individual lakes. I then calculated mean ratios 
for deep and shallow lakes and compared the means with a 
Mann-Whitney test. This method was used instead of analysis 
of covariance for three reasons: (1) significantly more 
variance was associated with the summer [TP] to spring [TP] 
relationship in Shallow lakes as compared to deep lakes 
(Elen t5rGPi<t07001),9(2) the data sets@wereasmall? andiXs) 
ratio data are often not normally distributed. Therefore 
parametric statistical techniques such as analysis of 
covariance, which assume homogeneity of variances and normal 


distributions, were inappropriate. 


C. Results and Discussion 

The relationship between spring [TP] and [Chla] was 
highly significant (P < 0.001) for both deep and shallow 
lakes (Fig. 1 and 2). The spring [TP]-summer [Chia] 


relationship for shallow lakes (Table 1) was: 
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logio{Chla] = 1.251 log, [TP]sp - 0.680 (r=0.83, n=25) (1) 


and for the deep lakes it was: 


log;o[Chla] = 1.015 log,o[TP]sp - 0.555 (r=0.80 n=31) (2) 


where [TP]sp is the [TP] at spring overturn. 
The correlation coefficients for these regressions (Eq. 


1 and 2) were not as high as for Dillon and Rigler's line: 


kogfiglGklali= if4£9ekog, ¢ITRlspe= Fii36 Mn=0n9Spen=46)ie (3). 


A partial explanation for why Eq. 1 and 2 had lower 
correlation coefficients than Eq. 3 is that the data sets 
for Eq. 1 and 2 covered a narrower range of spring TP 
concentraions than Dillon and Rigler's data set. Spring [TP] 
Values tcanged iirom Gisiadeto O152/ 4075. Dito 87959 2 and 7320 sto i180 
(nig/m*):sim Eg cbt a (2 ;'eand 13;2respectively: 

For analysis of covariance, I calculated a second set 
of regression lines. Only lakes with spring [TP] 
concentrations between 10 and 100 mg/m* were used. The 
spring [TP]-summer [Chla] relationship for the shallow lakes 


in this range was: 
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The equation for the lakes in this range from Dillon and 


Rigler (1974) was: 


logy, [Chiat =eiessoktogicITRispee, TeO052e(2=0390,"n=380) (3a). 


When Eq. la, 2a, and 3a were compared by analysis of 
covariance, the residual variances and slopes were not 
Significantly different (Table 2, Fig. 3). However, the 
adjusted means for Eq. la and 2a, and la and 3a were 
different, although the adjusted means for Eq. 2a and 3a 
were not different (Table 2, Fig 3.). Thus the spring 
[TP]-summer [Chla] relationships were different in deep and 
shallow lakes, but the two sets of deep lakes (Eq. 2a and 
3a) were Similar. Because of this difference in the spring 
[TP]-summer [Chla] relationship in deep and shallow lakes, 
Eq. 1 should be used to predict [Chla] in shallow lakes and 
Eq. 2 or 3 should be used for deep lakes. When more data are 
available, Eq. 1 should be compared with regression models 
for other shallow lakes to evaluate whether the spring 
[TP]-summer [Chla] relationship is consistent for shallow 
lakes. 

I did not attempt to incorporate total nitrogen (TN) 
(e.g., Smith 1982) into my models because my data set was 


too small for multivariate analyses. However, I had TN data 
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for one portion of my data set (Table 1, data from Prepas 
and Trew 1983). I used these data to test whether the mean 
summer [TN] to [TP] ratios were different in deep and 
shallow lakes. The ratios were different (P < 0.01). The 
mean [TN] to [TP] ratio was smaller in the shallow lakes 
than in the deep lakes (27 and 44, respectively) which 
according to Smith (1982) means the shallow lakes should 
have a lower summer [Chla] to spring [TP] ratio, than the 
deep lakes. Since the shallow lakes had a higher summer 
k@htaletouspringTiTPImravhor rhecencludedohIN<tokhnP] tratios 
in the shallow and deep lakes were not responsible for the 
differences found here. However, incorporation of [TN] into 
empirical models to predict summer [Chla] in deep and 
shallow lakes could reduce the residual variation. 

The higher elevation of the spring [TP]-summer [Chla] 
relationship for shallow lakes as compared to deep lakes 
(Fig. 3) supports the hypothesis that shallow lakes are more 
productive per unit of spring [TP] than deep lakes. The 
difference in elevation looks small when the lines are 
plotted together on a log-log scale. However when the 
regression lines are plotted in the arithmetic scale (Fig. 
4) the differences between the two models increase as spring 
[TP] increases. This trend suggests that shallow lakes 
respond more strongly than deep lakes to increases in spring 
nutrient concentrations. 

The differences between the predictions made by the 


spring [TP]-summer [Chla] models for deep and shallow lakes 
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have important implications for lake management. For 
instahcegrEq: leandi2 predicteas/Chiajtof' 28eand) 15 mg/m? , 
respectively for a lake with a spring [TP] of 50 mg/m°>. 
AccordingutoeDillon and! Riglery(4975);1a lake with a [Chla] 
of 15 mg/m° could support a warm-water fishery whereas a 
lake with a [Chla] of 28 mg/m* may not be able to. 

A partial explanation for the higher summer [Chla] to 
spring [TP] ratio in the shallow lakes as compared to deep 
lakes is the relationship between the spring [TP] and the 
average summer [TP] in the trophogenic zone of these two 
lake types. The mean ratios of summer [TP] to spring [TP] in 
Shallow and deep lakes were 1.57 and 0.87, respectively 
(Table 1). Thus [TP] increases in the shallow lakes and 
decreases in deep lakes from spring to summer. The 
difference between the mean summer [TP] to spring [TP] 
ratios was Significant (P < 0.001). Since phosphorus is 
generally the limiting nutrient in north temperate lakes 
(Schindler 1977) and [TP] increases in shallow lakes and 
decreases in deep lakes from spring to summer, the higher 
algal biomass per unit of spring [TP] in shallow lakes is 
reasonable. 

Clearly more information is needed on the cause of the 
increased summer [TP] in shallow lakes. Recycling of 
phosphorus due to lake mixing is one explanation for why the 
[TP] increases in shallow lakes (Oglesby and Schaffner 
1975). Phosphorus that would sediment out of the water 


column in a deep lake during the summer is resuspended in a 
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Shallow lake. 

Internal loading from the sediments is an additional 
explanation for the increase of [TP] in shallow lakes. Some 
authors suggest that internal phosphorus loading contributes 
Significantly to the phosphorus budget of shallow lakes 
during summer (Lie 1977; Larsen et al. 1981; Reynoldson and 
Hamilton 1982: Hanson and Stefan 1982; Landers 1982). 
However, the mechanism for internal loading to shallow lakes 
is still unclear. For example, Lie (1977) and Larsen et al. 
(1981) studied the same lake during the same time period; 
Larsen et al. reported that phosphorus release from the 
sediments during periods of temporary stratification 
accounted for most of the internal loading to the lake 
whereas Lie found phosphorus release from macrophytes was 
the major internal loading source. 

Clearly, more research on the processes controlling the 
[TP] in shallow lakes is needed. In particular phosphorus 
loading to shallow lakes needs better quantification. A 
clear understanding of these processes in shallow lakes will 
enhance our knowledge of these same processes in all lakes. 
It will also help to explain why the spring [TP]-Summer 


[Chla] relationship is different in deep and shallow lakes. 
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LPABLE Gs SDLINGdPOrewinter tObal. DhOoSpnorus concentration (LIP), 
summer chlorophyll a concentration. ([Chla]), summer [TP] in the 
trophogenic zone, and mixing regimes (D = deep, S = shallow) 

for the lakes used in the analyses. 


Lake Spring Summer Summer Mixing Reference 
[TP] [Chla] [TP] Regime 
(mg/m?) (mg/m?) (mg/m? ) 
Amnisk North 66.0 2ar6 3749 D Prepas 
Amnisk South 68.7 245. © 38.6 D and 
Baptist N. 57.56 34.1 65.6 S Trew 
Bapiass S. 713.9 50.9 554-3 D (1983) 
Birch 45.2 EE By; 5366 S 
Bourque 16.4 4.8 19.401 D 
Eden 28.7, 92 19.6 D 
Ethel 34.6 6.3 19.9 D 
Hastings 15226 36.6 156 S 
Hilda 28.4 5.6 25 tel S 
Hubbles 1980 27.4 10.4 Lo eo D 
Hubbles 1981 26.2 fons ZO ae D 
Marie $5.54 5.5 16.0 D 
Mink 24.7 Or 26.9 S 
Moore Sipe’ 4.7 23760 D 
Nakamun 1980 Soa Alege] 15 S 
Nakamun 1981 78.4 150 155 S 
Narrow 1981 tod 3.4 11.8 D 
Rol 1394 Sinise 14.8 D 
Sauer S540 13626 Cilhcow! D 
Star 2828 5.4 Zale 2 D 
Tucker 29 3030 58.9 S 
Twin 1980 21.4 7 eee) D 
Twin 1981 LG 23 126 D 
Wizard 3633 26.4 44,2 S 
Wolf 39.4 ipae 2423 D 
Lough Neagh Cae: 
1978 96 58 1 oe es S Gibson 
1979 98 100 jMBkay |e: S (pers. 
1980 114 79) 90 x S comm. ) 
1981 L234 Ose: 126. * S 
Oneida 
1975 174 PRA 87.8% S Boeri. 
1976 28.8 SPAS 25.3% S Mills 
1977 Sires 1228 eMOhase: S (pers. 
1978 262 6.9 Syiece S comm. ) 
1979 35.0 15 40.2 S 
1980 46.1 122.0 Givi e S 
198 1 44.2 TA<9 LOLA S 
1982 uae, 8.4 Bilvecet S 
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Lake Spring Summer Summer Mixing Reference 
[TP] [Chla] [TP] Regime 
(mg/m?) (mg/m?) (mg/m> ) 
Canadarago SOS DERE Tam Sa a S Oglesby 
Canandaiqua Nake a 2.6% ental D and 
Canadice O52 4.4+ sea D Schaffner 
Cayuga 1972 ZOvea Dhee ---- D (1978) 
Cayuga 1973 Di Ve Best siecle D 
Cayuga 1974 20i5 Sauee Sa D 
Conesus 1726 S. 65 tees D 
Hemlock LO BO Pers D 
Honeoye ees ire Zt i te S 
Keuka ‘Ebiea Sse SSeS D 
Lamoka ARS 272 Loot ae. 5 
Oneida 1973 42.0 29 S0rr aes S 
Otesego 8.4 et === S 
Otisco 8.4 Zee Se D 
Owasco 14.7 sh epiiaa Saleen! D 
Seneca eo 7. Sai ie D 
Skaneateles dom hope Shih ae D 
Waneta 24.0 ZO Or fe eked S 


* Integrated concentration for the whole water column. This 
data was not used to compare the relationship between summer 
PRP andespring [TP]. 


+ Summer [TP] of the trophogenic zone not available for 
Hubbles in 1980. Data for the [TP] of the epilimnetic zone was 
used instead. 


+ Chlorophyll a plus pheopigments. 
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Table 2. A comparison of three regression lines 
(Eq. la, 2a and 3a) by analysis of covariance 
(* indicates significant difference (P < 0.05)). 


COMPARISON OF THE THREE REGRESSION LINES 


Bartlett's Test for Homogeneity of Variances 


Rae 200 


Comparison of Slopes 


Pa 0.23 


Pairwise Comparison of Intercepts 
(Newman-Kuels test 


Shallow vs. Deep Pea O 0 Ones 
Shallow vs. Dillon and Rigler's P= 0,05) * 
Deep vs. Dillon and Rigler's P >> 01.150 
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Figure 1. Summer chlorophyll a concentration vs. spring 
total phosphorus concentration in the sample of shallow 
lakes. Dashed lines are the 95% confidence intervals for 
lakes outside the original data set. 
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Figure 2. Summer chlorophyll a concentration vs. spring 
total phosphorus concentration in the sample of deep 
lakes. Dashed lines are the 95% confidence intervals for 
lakes outside the original data set. 
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Eigurens. .oummegmchlorophyllwarconcentration vs. spring 
total phosphorus concentration for: (a) shallow lakes (Eq. 
la), (b) deep lakes (Eq. 2a), and (c) Dillon and Rigler's 
(1974) model (Eq. 3a) plotted together with the: (A) 


Original slopes, and (B) slopes averaged. 
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Figure 4. The regression equations of summer chlorophyll a 
concentration vs. spring total phosphorus concentration 
for the samples of deep and shallow lakes (Eq. 1 and 2) 

in the arithmetic scale. 
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III. The Role of Internal Phosphorus Loading in Two Shallow, 


Productive Lakes in Alberta, Canada 


A. Introduction 

During summer, total phosphorus (TP) concentrations 
increase in the surface water of many shallow lakes (Larsen 
et al. 1975; Stevens and Gibson 1977; Reynoldson and 
Hamilton 1982; Prepas and Trew 1983). These increases occur 
when external phosphorus (P) loading is low. Thus, in 
Shallow lakes, internal loading is often cited as the cause 
of the increased TP levels in the surface water. 

Three mechanisms for internal P loading to the surface 
water in shallow lakes have been suggested: (1) release of P 
from profundal sediments when the overlying water 1S anoxic, 
with subsequent transport of this P to the surface water 
(Larsen et al. 1981: Stefan and Hanson 1981); (2) 
resuspension of lake sediments during mixing (Reynoldson and 
Hamilton 1982); and (3) release of P from healthy and 
Scenescing macrophytes (Lie 1977; Carpenter 1980; Landers 
1982). 

Numerous laboratory studies have shown that P is 
released from lake sediments when the overlying water is 
anoxic (Mortimer 1941 and 1942: Theis and McCabe 1978; 
Holdren and Armstrong 1980) and release rates increase with 
higher water temperatures (Banoub 1975, Holdren and 
Armstrong 1980). Field studies have shown that P release 


from profundal sediments is an important source of loading 
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to the deep water of shallow lakes (Stevens and Gibson 1977; 
Larsen et al. 1981: Stefan and Hanson 1982). The link 
between deep water buildup of TP and fluctuations of [TP] in 
the surface water is less clear. Only one field study 
(Larsen et al. 1981) has shown that the [TP] increased in 
the surface water after mixing. However, in that study, TP 
increases in the surface water after mixing were 139-1000% 
greater than the losses from the deep water. To demonstrate 
that the P released from the sediments 1S transported to the 
Surface water during mixing, exchanges of TP between deep 
and surface waters need to be quantitatively balanced. 

Alberta has many shallow, productive lakes (Prepas 
1983). Total phosphorus concentrations in the surface water 
of these lakes increase during summer (Prepas and Trew 1983; 
Prepas and Wisheu unpublished data). To evaluate the 
contribution of the sediments to the TP levels in the 
surface water of these lakes, detailed studies were 
undertaken on Nakamun and Halfmoon Lakes. At Nakamun Lake, 
data were collected on the external TP loadings, vertical 
distributions of temperature, TP, total dissolved phosphorus 
(TDP), and dissolved oxygen (DO), and trophogenic 
concentrations of TP and chlorophyll a (Chla). At Halfmoon 
Lake, Similar data on temperature, TP, DO, and Chla were 
collected. Both lakes were sampled throughout the ice-free 
season. Also, sediment cores from Nakamun Lake were taken to 
the laboratory and incubated under oxic and anoxic 


conditions at three temperatures. These data were used to: 
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(1) evaluate the effect of water temperature and DO levels 
in the water overlying the sediments on P release rates in 
Shallow Alberta lakes; (2) quantify the relative importance 
of internal and external P loading to these lakes; (3) 
compare sediment release rates meaSured in situ and in 
cores; and (4) quantify the exchange of TP between the deep 
and surface waters during mixing and determine the effect of 


this input on TP and Chla levels in the surface water. 


B. Description of the Study Lakes 

The two study lakes, Nakamun and Halfmoon, are located 
in the boreal mixed-wood biome in central Alberta. Nakamun 
Lake is 70 km northwest of Edmonton and Halfmoon Lake is 20 
km east of Edmonton. Nakamun is a shallow (mean depth, 4.5 
m)(Fig. 1A), medium-sized lake (surface area, 3.5 km?), 
draining an area of 48.2 km’?, 42% of which is used for 
grazing and cereal crops. The mean depth of Halfmoon lake is 
Similar (4.8 m)(Fig. 1B), but the surface area is smaller 
(O42kme) peandethemdrainagesbasin,f57%a0fiwhich Ystused for 
grazing and cereal crops, is the same relative size (3.3 
km?). Both Nakamun and Halfmoon Lakes are eutrophic, the 
average summer Chla concentrations in 1981 were 150 and 46.2 
mg/m*?, respectively (Prepas and Trew 1983). Both lakes have 
relatively small littoral zones (< 20% of the sediment area) 
but the macrophyte beds in these areas are dense. Baseline 
chemical information on both lakes is in Prepas and Trew 
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C. Materials and Methods 

To determine how oxic and anoxic conditions in the 
water overlying the sediments and water temperature 
influenced P release from the sediments of Nakamun Lake, 
sediment cores were collected with a multiple corer 
(Hamilton et al. 1970). At the lake, all but the top 20 cm 
of sediment was released from the bottom of the core and the 
tubes were filled with lake water. The cores were 
transported to the laboratory and incubated at one of three 
EemMperatureseisee10, 25°C) under oxic or anoxic conditions 
for periods of 14 to 21 days. To maintain oxic (DO > 4 mg/L) 
and anoxic (DO < 1 mg/L) conditions, the water in the tubes 
was bubbled with air and nitrogen gas, respectively. 
Residual oxygen in the nitrogen was removed by passing the 
gas through a column of Catalyst R3-11 (Chemical Dynamics 
Corporation). To evaluate changes in TP levels in the water 
overlying the sediments, duplicate 5 mL subsamples were 
removed every three to five days and made up to 50 mL with 
double distilled water and stored for analysis. Water 
removed for analyses was not replaced. 

Bathymetric maps for both lakes were constructed from 
data collected with a Furuno model FE-400 depth sounder. 
Lake surface areas were determined from aerial photographs 
and lake watershed areas from 1:50,000 topographic maps and 
aerial photographs. Lake volumes were determined from 
bathymetric maps with a Tallos digitizer connected to a 


Hewlett-Packard model 9825B desk-top computer. 
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The TP loading to Nakamun Lake from terrestrial sources 
and losses via the outflow were determined from discharge 
and [TP] in the six inlet and the outlet streams. Stream 
discharge was determined by the dye dilution technique 
(Dillon 1974); the dyes used were Flourescein and Rhodamine 
WT. Water samples for TP analysis were collected midstream 
and placed in a 1-L Nalgene bottle. Loadings for the 20% of 
the watershed area with no streams, were estimated from 
runoff coefficients (kg/km*/day) from adjacent watershed 
areaS with streams. 

Aeolean P inputs to Nakamun Lake were estimated with a 
collector located in the lake, on a raft. The collector 
consisted of a Nalgene funnel leading into a 2-L Nalgene 
bottle (Gomolka 1975). To exclude insects, 250-um netting 
was secured to the bottom of the funnel (Dillon and Rigler 
1974b). To discourage birds from perching on the funnel, 
plastic straws were taped to the outside rim. Every one to 
two weeks the bottle was taken to the laboratory, the volume 
of water was recorded, and the water was analysed for TP. 
Aeolean samples collected in May and June were fouled by 
insects, birds, and humans. For this period, data collected 
for Wabamun Lake (25 km southwest of Nakamun Lake) were used 
to estimate aeolean loading (P. Mitchell, pers. comm.). 

Nakamun Lake was visited every three or four days and 
Halfmoon Lake was visited every two weeks during the 
ice-free season (between 11 and 13 May, and 10 and 15 


November, 1982). Vertical profiles of TP were determined on 
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water samples collected at 1i-m intervals from the surface to 
the sediment water at up to three stations in Nakamun Lake 
(8.5, 6.5 and 4.5 m deep; Fig. 1A) and over the deepest part 
of Halfmoon Lake (8.5 m deep; Fig. 1B). More than one 
Station was sampled at Nakamun Lake to evaluate if distance 
from the sediments was an important factor determining the 
TP level. At stations 1 and 3 in Nakamun Lake (Fig 1A), 
Water Samples fromedepths offd, 3) 5) 16¢tand Baand tprsyos, 
and 6 m, respectively were analyzed for TDP. At all three 
Stations in Nakamun Lake and at the one station in Halfmoon 
Lake, DO was determined on samples collected at i-m 
intervals from the sediment-water interface until 
well-oxygenated water (DO > 4 mg/L) was encountered. Water 
Samples for vertical profiles were collected with a 1.5-L 
drop-sleeve aluminum water bottle. Vertical profiles of 
temperature were meaSured at stations 1 and 3 in Nakamun 
Lake and at the one station in Halfmoon Lake at 1-m 
intervals with a Montedoro Whitney (model TC-5C) thermistor 
thermometer accurate to 0.1 °C. The epilimnion was defined 
according to Hutchinson (1957). The trophogenic zone was 
estimated as two times Secchi disk depth (Dillon and Rigler 
1974a), Trophogenic TP and Chla levels representative of the 
whole lake were collected at 5 to 10 stations in Nakamun and 
Halfmoon Lakes. Samples were collected from the surface to 
the bottom of the trophogenic zone with a polyethylene tube 
and pooled in 2-L opaque Nalgene bottles. In Nakamun Lake, 


this water was also used to determine the phytoplankton 
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Species composition. All Nalgene bottles used for 
transporting samples back to the laboratory were treated as 
outlined in Prepas and Rigler (1982). 

Total phosphorus was determined on 50 mL subsamples 
with the potassium persulfate technique described by Prepas 
and Rigler (1982), with one modification. Since the 
Aphanizomenon colonies in Halfmoon Lake did not pass through 
a 250-um net, the samples were not filtered. Water for TDP 
analysis was filtered through a pre-rinsed 0.45-um Millipore 
HAWP membrane filter and analysed as for TP. Duplicate 
Subsamples were analyzed from water over the sediment cores 
and for TDP; the remaining analyses were in triplicate. 

Dissolved oxygen was determined on water samples which 
were fixed in the field and analyzed by Carpenter's (1965) 
modified Winkler technique. For the water overlying the 
sediment cores, DO was determined by a micro-Winkler 
technique (Burke 1962). Chlorophyll a was determined with 
the ethanol extraction method described in Bergmann and 
Peters (1980), with one exception. I used a Bausch and Lomb 
Spectronic (Spec) 710 spectrophotometer instead of a Spec 
100. When Chla samples were read on both machines, the Spec 
710 consistently gave higher readings. To make my data 
comparable with those reported from Spec 100s, my Chla 
values were corrected based on a comparison using Chla from 


Sigma Chemical Co.: 
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[Chla](Spec 100) = 0.88 ([Chla](Spec 710)) + 0.28 
Un = 20.965) neeee P< 10.01) 


Phytoplankton counts were determined by the inverted 
microscope method described by Lund et al. (1958). 

To estimate stratum and whole-lake masses of TP, 
profiles were collected at three stations in Nakamun Lake. 
The lake was divided into three sections corresponding to 
the depths of Stations 2, 3, and 1: areas over maximum 
Gdepensmors0 to 5 Mio to / im, andy / £0 8.5m, respectively. 
Masses of TP were calculated for each Stratum, over each 
section, and summed for whole-lake masses. Masses of TDP 
were calculated in a similar manner, except the lakes was 
divided into two sections: areas over maximum depths of 0 to 
Deans eb O loo om. 

To calculate P loading to the deep water from the 
sediments and to the surface from the deep waters, the lakes 
were divided into deep and surface strata. The division was 
at 2.5 m because water above this depth was always well 
oxygenated and the trophogenic zone was generally in this 
region (average depths of the trophogenic zone during the 
Summer for Nakamun and Halfmoon Lakes were 1.6 and 2.0 m, 
respectively). 

The net external loading (NEL) to Nakamun Lake between 
two sampling dates, t to ttn (where mM is the number of 


days), was calculated with the mass balance equation: 
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NEL = TPin - TPout (iy) 


where TPjn is the TP entering the lake from runoff and 
aeolean loading, and TPoyt is the amount of TP lost from the 
Wake viarthe=outriowlfrom’izto f+m. “External loading to 
Halfmoon Lake was estimated from coefficients calculated 
with data from Nakamun Lake. Net internal load (NIL) to 
Nakamun Lake from ft to ft+n was calculated from the mass 


balance equation: 


NIL = /\TP - NEL (2) 
where /\TP is the change of TP in the whole lake from ft to 
ttn. The units for NEL and NIL are kg. During periods when 
the lakes were stratified, release rates of TP and TDP (RRTp) 
from the sediments (mg/m?/day) were calculated with the 


equation: 


RRPp*EeTPraey) /* (ne/ vas) (3) 


wheres?Pre)usyene Ehange m@TPeor*TDPM(ng)ebelow 275%m'/9n 
is the number of days the lake was stratified, and As is the 
surface area (m?) of the sediments which were anoxic. 
Release rates from the cores were calculated with a modified 
EQumo tet DwaAvewasetheeinereace OLelTP=an the watereabove ‘the 
core, NM was the number days this water was bubbled with 


nitrogen gas, and As was the area inside the core (11 cm’). 
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Thermal budgets were constructed for Nakamun Lake to 
provide independent estimates of TP movement from the deep 
water to the surface water during mixing events. The volume 
of the epilimnetic water required to raise the temperature 
of a deep water stratum (e.g., 4.5 to 5.5 m) from that 
measured at t to that measured at (ft+n) (V_.;.) was 


mix 


calculated: 


if a 
Vote = Vd ( Fote stat oun Td ) (4) 


opie cdit 
where Vy is the volume (m?) of the deep water stratum, Td,t 
and Ty +4, are the temperature (°C) of the stratum on t and 


t+n, respectively, and T is the volume-weighted 


epi 
temperature of the epilimnion at ft. The estimated volume of 
water exchanged between t and t+n and the TP concentrations 
at time tf in the two strata, were used to estimate the 
vertical flux of TP between ft and ttn. This flux was 
estimated for each stratum below the eplimnion where the 
temperature had increased from t to ft+n. The total flux was 
used to calculate a predicted [TP] in the surface water at 
ttn. These budgets assume that during mixing, all exchange 
was between the deep water and the epilimnion (i.e., water 
did not mix between the deeper strata). Thermal budgets were 
calculated for periods of mixing in which the thermal mass 


(temperature x volume) of the lake did not change 


Significantly (< 5%) between t and t+n. 
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D. Results 


Sediment Core Experiments 

In the sediment core experiments, TP was released when 
the overlying water was anoxic (Fig. 2). However, when the 
Same water was aerated, TP levels decreased in the water 
column (Fig. 2). Release rates (RRyp) increased proportional 


to water temperature (Table 1, Fig 2B): 
RR = 0.51 @ + 1.78 (r = 0.98, n = 11, P < 0.001) (5) 


where @ is temperature (°C). These results suggest that: (1) 
when the water overlying the sediments in Nakamun Lake is 
anoxic, TP will increase in this water; and (2) in 
individual lakes, release rates are temperature ea nay 
Although core experiments were not run on sediments from 
Halfmoon Lake, I expected similar results since both lakes 
have similar mean depths (Fig. 1), superficial sediment type 
(black, organic material), levels of productivity, and water 


chemistry (Prepas and Trew 1983). 


Phosphorus Loading to Nakamun Lake 

Nakamun Lake had 12 successive periods of water 
stagnation and mixing between 17 May and 2 November, 1982. 
These periods of stagnation and mixing were divided into 
three categories: (1) stratification, when the water next to 


the sediments was anoxic and the chemocline or thermocline 
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was not depressed from the previous sampling day; (2) 
complete mixing, when chemical or thermal stratification was 
not evident; and (3) partial mixing, when the thermocline or 
chemocline was depressed from the previous sampling day, but 
the water over the sediments was anoxic. 

Nakamun Lake was stratified four times between 25 May 
and 23 August (Table 2); one of these periods lasted 34 days 
and the other three lasted 6 or 7 days. During these 
Stratified periods, DO was rapidly depleted in the water 
over the sediments (Fig. 4B). For instance, on 22 July, the 
[DO] at 8 m at Station lewas 5:4 mg/L; on 26 July, the [DO] 
was 0.3 mg/L at the same location. As DO levels approached 0 
mg/L, the TP and TDP concentrations increased near the 
sediments (Fig. 4C and D), TP and TDP mass increased in the 
deep water, and the net internal load (Eq. 2) was positive 
(Table 2, Fig. 3A). Total phosporus levels in the surface 
water did not change more than 5% during any stratified 
period (Table 2). 

During seven periods from 17 May to 2 November, Nakamun 
Lake mixed completely or partially. To facilitate the 
analysis, the period from 23 August to 2 November was 
divided into two periods (23-26 August and 26 August-2 
November). Thus, Table 2 lists eight periods of complete or 
partial mixing. During the two periods of complete mixing 
which directly followed other periods of complete mixing 
(17-25 May and 26 August-2 November) (Table 2), the [TP] was 


constant throughout the water column (Fig 4C), the water 
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over the sediments was well oxygenated (Fig. 4B), and the 
net internal load was negative (Table 2). During these two 
periods, the loss of TP from the whole lake was 6 and 
2%/day, respectively. The remaining six periods of complete 
Or partial mixing were between 28 June and 26 August (Table 
2). Prior to each of these six periods, the water next to 
the sediments was anoxic and TP and TDP levels were high in 
the deep water (Table 2, Fig. 4). When the lake mixed, the 
P-rich deep water exchanged with the relatively dilute 
Surface water. Consequently, during five of these periods, 
TP increased 3-43% in the surface water (Table 2, Fig 3A). 
During the sixth period (9-16 August), the TP level in the 
surface water decreased. Total dissolved P levels in the 
surface water remained constant throughout the summer (Fig 
4D). During mixing events, TDP which was transported to the 
surface water, was probably incorporated immediately into 
the particulate P pool by the phytoplankton (Lehman and 
Sandgren 1982). 

There were four periods of complete mixing between 15 
July and 26 August (Table 2). During three of these periods, 
both TP and TDP levels decreased in the deep water and the 
net internal load was negative (Table 2). During these 
periods, the loss of TP from the deep water was greater than 
the increase in the surface water (Table 2). Thus, some of 
the P lost from the deep water went to the sediments. During 
the fourth period of complete mixing (23-26 August) the TP 


increased in both the deep and surface water and the net 
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internal load was positive (Table 2). During the two periods 
of partial mixing in Nakamun Lake, TP and TDP continued to 
increase in the water over the sediments (Fig. 4C and D). 
Consequently, losses from the deep water were less than 
increases of TP in the surface water during the first period 
of partial mixing (28 June-5 July) and TP increased rather 
than decreased in the deep water during the second period of 
partial mixing (Table 2). During these two periods, the net 
internal load was poSitive. 

Terrestrial loading to Nakamun was low, 39 kg between 
17 May and 2 November. When the terrestrial loading 
estimates collected between 17 May and 2 November were 
extrapolated to a full year, the coefficient for the 
watershed was 1.9 mg/m?. Between 17 May and 2 November, 
aeolean sources contributed a similar amount of TP to the 
lake, 43 kg. When the lake was stratified, the net external 
load (Eq. 1) was insignificant compared to the net internal 
load (Eq. 2) (37 and 1468 kg, respectively). During the four 
periods of stratification and the second period of partial 
mixing (5-15 July), release rates of TP and TDP from the 
sedaments rangedifromei0TS tos 19<6j andr 8s if tom28.0 


mg/m*/day, respectively. 


Phosphorus Loading to Halfmoon Lake 
Halfmoon Lake was thermally stratified in late May and 
by early July, partial mixing was evident (Fig. 5A). 


However, the mixing in Halfmoon was weaker than in Nakamun 
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Lake. Halfmoon Lake did not completely mix until 25 
September. Consequently, the water overlying the sediments 
in Halfmoon Lake was anoxic from 27 May until 25 September. 
Based on thermal patterns, the study on Halfmoon Lake was 
broken into five periods: (1) a stratified period; (2) and 
(3) two consecutive periods of of partial mixing (the first 
was much weaker); and (4) and (5) two consecutive periods of 
complete mixing (the first followed a TP buildup in the deep 
water). 

During the period when Halfmoon Lake was stratified (27 
May-25 June), TP levels decreased 2%/day in the surface 
water and TP increased (78) kg in the deep water (Table 2, 
Figt3Bis Duringhthe.firstnperiodwof partial imixingn (25 
June-8 August), TP increased 36% in the surface water as the 
result of P-rich deep water mixing with the relatively 
dilute surface waters (Table 2, Fig. 3B). During this 
period, TP continued to increase in the deep water since the 
water overlying the sediments remained anoxic (Table 2, Fig. 
3B, 4B and C). Mixing was stronger from 8 August to 10 
September. Total phosphorus increased 52% in the surface 
water and a net loss of TP was registered Fon the deep 
water (Table 2, Fig 3B). However, the mixing was still 
incomplete during this period and TP levels remained high in 
the deep water (Fig 5C). From 10 to 25 September, the lake 
mixed completely, the TP resevoir in the deep water was 
distributed throughout the water column, the TP level in the 


Surface water increased 31%, and the TP level of the deep 
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water decreased (Table 2, Fig 3B and 5). During this period, 
TP increases in the surface water accounted for only one 
third of the losses from the deep water (Table 2). The lake 
continued to mix from 25 September to the end of the study 
(29 October). Total P decreased during this period, at a 
rate of 1%/day (Table 2). 

When the water next to the sediments in Halfmoon Lake 
was anoxic, the net increase of TP was 147 kg (Table 2). 
Runoff and aeolean loading coefficients for Nakamun Lake 
(0.0051 and 0.065 mg/m*/day, respectively) were used to 
estimate external loading to Halfmoon Lake. From 27 May to 
10 September, external loading supplied an estimated 5 kg of 
TP to Halfmoon Lake. AS in Nakamun, external loading was 
insignificant in Halfmoon Lake compared with internal 


loading. 


Thermal Budgets 

Thermal budgets were used to predict TP concentraions 
in the water above 2.5 m in Nakamun Lake on three occasions: 
1, 5, and 15 July. The predicted TP concentrations were all 
higher than the measured values (Table 3). However, the 
predicted TP values assume no sedimentation. When a 
sedimentation rate of 1.5%/day (Rigler 1974) was applied, 
the predicted and measured values were much closer (Table 
3). One predicted value was still overestimated (15 July) 
and another was underestimated (5 July). These calculations 


assume that the lake mixed immediately after the previous 
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sampling day. Thus, the calculations are only approximate, 
but they suggest that sedimentation rates were variable in 


Nakamun Lake. 


Chlorophyll a and Total Phosphorus in the Trophogenic Zone 
In both lakes, TP and Chla concentrations in the 
trophogenic zone were correlated on a day to day basis from 
May to: Augustes(r~= > 0 805! ns =.:28),0 P< 0.001c and m= 0.76, n = 

7, P< 0.05 for Nakamun and Halfmoon Lakes, respectively) 
(Fig. 6). As well in Nakamun Lake, TP and Chla levels were 
correlated during the fall (September to November, pr =0.98, 
n= 10, P = 0.001). However in Halfmoon Lake, TP increased 
in the fall due to lake mixing, while Chla levels decreased 
(Fig. 6). Thus, P transported to the surface water during 
the growing season (May-August) was incorporated by the 
phytoplankton, whereas P transported to the surface water 


during the fall was not. 


E. Discussion 

The data from Nakamun and Halfmoon Lakes support the 
Original hypothesis that P is released from lake sediments 
when the overlying water iS anoxic. During these periods, 
internal loading contributed 1468 and 147 kg of TP to 
Nakamun and Halfmoon Lakes, respectively. The sediments were 
likely the source of this internal loading: (1) the buildup 
of TP and TDP began and was always greatest at the 


Sediment-water interfaces(Fig. 4C, 4D, and 5C); (2) in 
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Nakamun Lake, where both TP and TDP were measured, the 
build-up in the deep water was TDP not particulate P (during 
the four periods of stratification, TDP were 77-226% of the 
increased TP in deep water); (3) in four out of five cases 
when the lakes were stratified, increases of TP in the deep 
water could not be accounted for by losses of TP from the 
Surface water (Table 2). Although groundwater is a potential 
P source, it was not monitered in either lake. However, 
loading patterns to the lakes indicate that the sediments, 
rather than groundwater, were the major TP source. If 
groundwater had been the major P input, the net internal 
load should have been positive throughout the summer, not 
just during stratified periods (Table 2). 

When Nakamun and Halfmoon Lakes were stratified, 
external loading was insignificant (37 and 5 kg, 
respectively) compared with internal loading (1468 and 147 
kg, respectively). These differences were not due to an 
underestimation of external loading. Aeolean loading 
coefficients from this study were similar to those for 
Wabamun Lake (25 km southeast of Nakamun) (P. Mitchell, 
pers. comm.) and Lake St. Nora, Ontario (Gomolka 1975). 


Aeolean TP Loading 


(mg/m? /day ) 
Nakamun Wabamun St. Nora 
12sJvle12"Aug 3.24 3. 08 2578 


23 Aug-23 Sep 1.43 2.20 24.03 
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The yearly terrestrial loading coefficient calculated 
for Nakamun Lake was lower than those measured for 10 other 
central Albertan (P. Mitchell, per. comm.) and 31 Ontario 
(Dillon and Kirchner 1975) watersheds with similar land use 
(1.9, 12.5, and 28.8 mg/m?/yr, respectively). Our 
coefficient was low, partly because spring runoff can 
account for up to 70% of terrestrial loading to local lakes 
(Mitchell and Hamilton 1982) and only the last part of 
Spring runoff was measured in my study. But even with a 
correction for spring runoff, my loading value was only half 
the amount for the other Albertan watersheds. However, this 
difference was too small to affect the relative difference 
between external and internal loading in Nakamun and 
Halfmoon Lakes. 

The in-situ TP release rates estimated for Nakamun and 
Halfmoon Lakes were higher than those predicted from the 
sediments incubated in the laboratory (Eq. 5) (Table 1). 
However, there are inevitable errors associated with 
extrapolating release rates estimated from a 11 cm? 
microcosm to whole-lake release rates: (1) there was no 
sedimentation of organic material to the water overlying the 
cores; (2) adsorption of P by biota on the walls of the tube 
housing the core (Rigler 1956); and (3) imprecise 
measurement of the sediment area of the lake. Considering 
the possible errors, the predicted and measured rates are 
reasonably close (Table 1). In-situ release rates of TDP in 


Nakamun Lake were less than for TP on three occasions and 
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more than for TP on two ocassions (Table 1). These 
differences are less than 40%, except for the period of 22 
to 29 July (Table 1). During this period the TDP release 
rates were 228% higher than the TP release rates. Small 
differences are expected since TDP wasS not measured at as 
many depths or stations as TP. However, the differences for 
the period of 22 to 29 July were due to a large loss of 
particulate P throughout the water column. 

Losses of TP from the surface water to the deep water 
were quite low in Nakamun Lake during the stratified periods 
(< 1%/day). The sedimentation rate was much higher in the 
mixing period of 17 to 25 May (6%/day) and somewhat higher 
in the mixing period of 26 August to 2 November (2%/day). 
One explanation for low sedimentation rates from the surface 
water during stratified periods, is that P can be 
transported from the deep, to the surface waters by eddy 
diffusion when vertical TP gradients are high (Larsen et al. 
1981). Lower sedimentation rates would also be expected from 
June to November than from 17-25 May in Nakamun lake because 
of differences in the algal community composition. Diatoms 
were the dominant algal group in May, when sedimentation was 
high. Whereas cyanophytes, which have lower sedimentation 
rates, were dominant for the rest of the study (Fig. 7) 
(Reynolds et al. 1982). 

This study supports the hypothesis that much of the P 
which is released from the sediments under anoxic conditions 


is transported to, and remains in, the surface water when 


25 36 bswwiaem ton, Sew. Sar es 
101 eentezstiiy — 
to) eeol sprit: ov subs 
sneules rer 

1936¥ qeob ort - tetew a2n 
ebotisag beltigente aie. ant305 Suk sumevan oH yo a ati W 
sit od sotpdnd aioum ene end iteaioniae tae i) 
redgid texwesioa Brie (vad \ta) ca ae ed 1 Bebodene pe 
.(yab\yaS) aeditevoK ¢ of Saupua dS to- foleserntese a) 


eosiwe od5 wet) ¢8967 sotretmemibee wot: sot oe a at 


et pds Sari ie painter cont 
ghbs yd. eupsew onelcbe. oda) Gd vetoeb a 


{8 36 Aeated) deid sxe etastherp 9 Ree sega | 
moi betseqes sad. odie bicew! gepag pemmpbhel a 
seysred stad numeiet af ene de+Tr trond alienate ee 

amossid: .nolpzigogmes, asian, Sin oe | ib 38 Di 
saw noldscoamibee tiedy | ysM ob weg mee stow 
acitesnemibse zevoL avad stoldw oi 2) Bi 

(« .pi9) Aibuate, ots so see a 0 | % pres x 


: =) mie ae. 
@ sf2) do de zi perigee 


ce 


48 


the lake mixes. In both lakes, TP increased in the surface 
water during eight of nine mixing events, which followed a 
build-up of TP in the deep water. During periods of complete 
mixing, increases in the surface water could not be 
accounted for by external loading. However, these increases 
could be accounted for by a transfer of TP from the deep 
water, in all but one case. During periods of partial 
mixing, increases in the surface water were greater than 
losses from the deep water because P release from the 
sediments continued. Independent estimates of TP movement 
between the deep and surface waters (thermal budgets) also 
indicated that lake mixing caused the TP increases in the 
Surface water. For most of the study period, macrophyte 
release and resuspension of sediments during mixing do not 
appear to be important P sources. If these sources were 
important, net internal P loading should have been positive 
rather than negative during periods of complete mixing. 
However, during one period of complete mixing in Nakamun 
Lake, the net internal load was positive rather than 
negative (23-26 August) (Table 2). During this period, TP 
increased in the surface water (as expected during mixing), 
but a corresponding decrease was not registered in the deep 
water. The P input to Nakamun lake for this period could 
have come from a source other than sediment release, such as 
senescence of macrophytes in the littoral zone (Carpenter 


1980; Landers 1982). 
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The data from Nakamun and Halfmoon Lakes clearly 
illustrate the importance of P release from the sediments on 
the TP budget and [TP] in the surface water in shallow 
lakes. The data also illustrates the importance of the 
increased TP levels in the surface water on algal biomass 
(Fig. 6). Further studies are needed on the factors that 
control TP and Chla levels in shallow lakes on a year to 
year basis. For example, summer Chla concentrations in 
Nakamun Lake in 1980 through 1982 were 42, 150 (Prepas and 
Trew 1983), and 60 mg/m’, respectively; and for Lough Neagh, 
Ireland in 1978 through 1981 average summer Chla 
concentrations were 58, 100, 77, and 103 mg/m’, respectively 
(C.E. Gibson pers. comm.). These differences may be due to 
variations in mixing patterns and release rates from year to 
year. For instance, had strong mixing occurred during June 
instead of July in Nakamun Lake, temperatures at the 
sediment-water interface would have been higher and more P 
would have released into the deep water. When the lake 
mixed, TP increases in the surface water would have been 
even greater and Chla concentrations may have been higher. 
The importance of DO levels and temperature on P release 
rates were also demonstrated in this study. Further studies 
are needed on the factors which control between lake 
differences in P release rates. As the factors controlling 
year to year variation in productivity and lake to lake 
variation in P release rates are established, management 


tactics for shallow lakes will improve. 
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Table 1. A comparison of release rates of total phosphorus 

(TP) from Nakamun Lake sediments incubated at three 

temperatures in the laboratory and from in-lake TP budgets 

for Nakamun and Halfmoon Lakes during periods of chemical 
Stratification in 1982. Measured (Meas) release rates of TP 

for both lakes are compared with release rates predicted from 
incubated sediment cores (Eq. 5). The difference between 
measured and predicted (Diff) and the measured release rates of 
total dissolved phosphorus (TDP) for Nakamun Lake are also given. 


Temperature Release Rate (mg/m?/day) 


ee TBP 

Date oC Meas Pred Diff Meas 
Mart Cores 3 SiG 
Apr Cores 10 Ge 
Mart Cores 25 14.6 
25 May-28 Jun Nakamun 14 1 O19 Og tes 6 ot 
OSedul=toe duly iNakamun V7, Lao 1Ore wits | ate 
22eJ001-29 sully sNakamun 17 tes Me segs) Lege) 20U 
03 Aug-09 Aug Nakamun 19 1407 Pleo Sipe 8.9 
16 Aug-23 Aug Nakamun 1 AZo MOREE pel 16.0 
27 May-25 Jun Halfmoon 10 18e2 6.9 625 
11 Jun-08 Aug Halfmoon ii 1 epee Tee oe 


+ These cores collected in 1983. 
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Table 3. A comparison of measured and 
predicted total phosphorus concentrations 
((TP]) imsthée. top, 2.5¢m of-the-water. column of 
Nakamun Lake after periods of mixing. 

Predicted values were calculated from thermal 
budgets (Eq. 4) assuming no sedimentation from 
the zone above 2.5 m (0% Sed) and sedimentation 
Fates O07 1.5%/day G1. 5% Sed):. 


[TP] (mg/m?) Above 2.5 m 


Measured Predicted 
Date 0% Sed “175% "Sed 
Otsu 45 80 75 
05 dud 94 97 9 1 
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Plgure |. Bathymetric maps of: (A). Nakamun and (B) Halfmoon 


Lakes. 
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TOTAL PHOSPHORUS (mg/m?) 


TIME (DAYS) 


Figure 2. Mean total phosphorus concentration (+ SE) of the 
water over lake sediment cores versus time. Cores were 
incubated in the dark under: (A) anoxic and oxic conditions 
at 10° C and (B) under anoxic then: oxic conditions 
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Figure 3. Mass of total phosphorus above (>) and below (<) 
the depth of 2.5 m versus time in: (A) Nakamun and (B) 
Halfmoon Lakes. Mass of total dissolved phosphorus (TDP) 
below 2.5 m also given for Nakamun Lake. Black bars indicate 


periods, Of Stratification. 
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Figure 4. Three dimensional graphs of time versus depth 
versus: (A) temperature, (B) dissolved oxygen concentration 
({pD0O]), (¢€) total» phosphorus concentration ({TP]), and 

(D) total dissolved phosphorus concentration ([TDP])over 
the deepest station at Nakamun Lake. Note that the depth 
axes are reversed in C and D. 
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Figure 5. Three dimensional graphs of time versus depth 
versus: (A) temperature, (B) dissolved oxygen concentration 
({DO])), and (C) total phosphorus concentration ([TP]) 

in Halfmoon Lake. Note that the depth axis is reversed in C. 
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[TP] and [Chla] (mg/m?) 
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[TP] and [Chia] (mg/m?) 


Figure 6. Total phosphorus concentration ([TP]) and 
chlorophyll a concentration ([Chla]) in the trophogenic 
zone from May until October, 1982 in: (A) Nakamun 

and (B) Halfmoon Lakes. 
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Fagurery .} Percents composition of,dominant .groups.in the 
algal community in the trophogenic zone of Nakamun Lake 
from Mayountodl October: h982: 
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IV. Concluding Discussion 

A comparison was made of the relationship between 
Spring total phosporus (TP) and summer chlorophyll a 
concentrations in deep and shallow lakes. The relationships 
were different (P < 0.05) in the two lake types. Shallow 
lakes had a higher summer [Chla] per unit of spring [TP] 
than deep lakes. The suggested reason for this difference 
was that TP concentrations in the surface water increased 
from spring to summer in shallow lakes, but decreased in 
deep lakes. Thus, relatively more phosphorus (P) was 
available to the phytoplankton community in shallow lakes 
during the summer than was indicated by spring [TP]. 

Internal P loading from the sediments is often 
suggested as the reason for increased surface water 
concentrations of TP in shallow lakes. The importance of 
this source was investigated in two shallow lakes in 
Alberta, Nakamun and Halfmoon. Laboratory experiments on 
sediment cores indicated P released into the overlying water 
when this water was anoxic. When the lakes were stratified, 
dissolved oxygen levels in the water overlying the sediments 
were low and TP levels increased in this water. During these 
periods, internal loading contributed 1468 and 147 kg of TP 
to Nakamun and Halfmoon Lakes, respectively, while external 
loading only contributed 37 and 5 kg, respectively. Total 
phosporus release rates, calculated from in-lake budgets for 
Nakamun and Halfmoon (12.7 and 15.6 mg/m?/day, 


respectively), were higher than those predicted from cores 
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incubated in laboratory (9.7 and 7.2 mg/m?/day, 
respectively). Considering the possible errors in 
extrapolating laboratory measurments to the field, the 
predicted and meaSured rates were reasonably close. Lake 
mixing transported the P which released from the sediments 
to the surface water. After eight of the nine mixing events 
which immediately followed stratified periods, the [TP] 
increased 3-43% and 31-52% in the surface water of Nakamun 
and Halfmoon Lakes, respectively. Independent estimates of P 
movement in Nakamun Lake (thermal budgets) confirmed that TP 
increases in the surface water were due to lake mixing. The 
Subsequent increase of TP levels in the surface water 
resulted in increased Chla levels. 

Although sediment release was the major source of 
internal loading, there were other internal inputs of P in 
Nakamun Lake. The importance of other mechanisms of internal 
P loading to shallow lakes (macrophytes (Lie 1977; Landers 
1982) and sediment resuspension (Reynoldson and Hamilton 
1982)) should be studied in conjunction with release from 
the sediments. Further studies are also needed on internal P 
loading from the sediments in other shallow lakes. From 
these studies, the effect of different mixing patterns on P 
release and transport can be determined. For instance, will 
more P be transported to the surface water in lakes which 
mix completely during the summer as compared to lakes which 
only partially mix? More studies on internal loading from 


sediments will also determine if P release rates are 
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Significantly different between lakes. Since internal 
loading seems to be the important factor influencing summer 
TP levels in shallow lakes, more productive lakes may have 
higher release rates. As the relative importance of mixing 
patterns and release rates are eStablished, management 


tactics for shallow lakes will improve. 
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V. Appendix A: A Test of Two Chlorophyll Models With An 
Independent Data Set. 

Empirical models to predict the growing season mean 
concentration of chlorophyll a ([Chla]su) are based on 
phosphorus (P) (Nicholls and Dillon 1978). The relationship 
between total P concentration ([TP]) and ([Chla]su) is quite 
robust (Prepas and Trew 1983). However, the confidence 
limits for the models are quite broad (Dillon and Rigler 
1974; Nicholls and Dillon 1978). Smith (1982) attempted to 
improve the predictability of [Chla]su by developing a model 
based on two variables, TP and total nitrogen (TN) 


concentrations: 


log, ob¢hiaisu =m0t653 LogpalTRisuhhs0 2548 thogrebTNisulahiedS17 (1) 


where [TP]su and [TN]su are mean growing season [TP] and 
[TN], respectively. The model is based on 311 cases. 

In order to demonstrate how much variability TN 
accounts for in the phosphorus-chlorophyll relation, Smith 
compared his model (Eq. 1) with Dillon and Rigler's (1974) 


modelsabased ome[TP}talone: 


lege. hehiadisur-or /44¢9cho0q%>5 LEP] Sp? jnp1in36 (2) 


where [TP]sp is the lake's [TP] at spring overturn. For the 


test, Eq. 1 and 2 were used to predict [Chla]su for 311 
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lakes based on [TP]su and [TN]Jsu (Smith 1982). Thus, [TP]su 
was substituted for [TP]sp in Eq. 2. To compare the accuracy 
of the models, residual variations (RV) were calculated for 


each model: 


RVie KGobsumseCoredia/dne zed (3) 


where Cobs is the observed [Chla]su and Cpred is the 
predicted [Chla]su. For the 311 cases, the RV for Eq. 1 was 
1% of the RV for Eq. 2. When only those lakes with a TN to 
TP ratio > 12 (i.e., lakes which are phosphorus limited 
according to Dillon and Rigler (1974)) the RV for Eq. 1 was 
6027 of ‘thet RVeior  Eqardr. 

However, Smith's comparison has a major problem (which 
he recognized in his paper). Equations 1 and 2 were compared 
with the same 311 cases used to construct Eq. 1. Thus, it 
was inevitable that Eq. 1 would be a more accurate predictor 
of [Chla]su than Eq. 2. The two models should have been 
tested with an independent data set. 

To evaluate both models, I used an independent data set 
from western Canada consisting of [TP]sp, [TP]su, and [TN]su 
(Prepas and Trew 1983). For the comparison, I used 20 of the 
22 lakes which had sufficient data for this comparison. Two 
lakes were excluded @ priori because the [TP]sp in these 
lakes (303 and 761 mg/m*) were well beyond the range in 
Dillon and Rigler's study (3 to 180 mg/m*). Three sets of 


predictions were made: (1) Both [TP]su and [TN]su were used 
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in Eq. 1 to predict [Chla]su, (2) [TP]su was used in Eq. 2 
to predict [Chla]su, and (3) [TP]sp was used in Eq. 2 to 
predict [Chla]su. Contrary to Smith's approach, I calculated 
RV's for each model based on standard statistical 
techniques. Thus, the denominator in Eq. 3 was replaced with 
n- i where i is the number of parameters estimated in the 
corresponding regression model (i.e., / = 3 and 2 for 
predictions from Eq. 1 and 2, respecitively). 

The RV's from the three sets of predictions were 
compared pairwise with an F-test. Since three comparisons 
were made, the chance of committing a type I error was 
greater than the value given in a standard F table (Zar 
1974). Therefore, a probability level of 0.02, rather than 
the traditional 0.05, was used to distinguish significance. 

For the western Canadian lakes, Eq. 2 with [TP]su as 
the independent variable was the best predictor of [Chla]su, 
followed by Eq.'s 1 and 2 where [TP]sp was the independent 


Variable (Table 1). The accuracy of Eq. 1 and 2 was 
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Statistically indistinguishable. (F Voie P= 
Oho utcr Bo. 2 basedmon 2P lisp; — = 12.65, 0f- =" 17,16, 2 = 
0.04 for Eq. 2 based on [TP]su). Thus for the 20 lakes from 
western Canada, the addition of [TN]su does not improve the 
predictability of [Chla]su over [TP]su or [TP]sp alone. 
However, Eq. 2 was a better predictor of [Chla] when [TP]su 
WAS) USeC rather thane tupien (fF = 94-8, 7dr = 18, 18, P =0.002).. 


These results show [TP]su is a better predictor of [Chla]su 


than [ETP Isp, 
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There are two possible reasons why the addition of TN 
did not add significant information to the [TP]-[Chla]su 
relationship for the 20 lakes in western Canada: (1) All 
these lakes had TNsu to TPsu ratios > 18 (Prepas and Trew 
1983). Consequently, nitrogen was probably not limiting 
algal production (Dillon and Rigler 1974), although Smith 
argues that nitrogen is limiting even when the TNsu to TPsu 
ratio is > 18. The*addition’of nitrogen*to the [TP]-[Chla]su 
model may be significant in lakes where the TNSu to TPsu 
ratio is < 18 and should be tested with the appropriate data 
set. (2) The data set 1S quite small and only from one 
region. Equation 1 should be tested on a larger data set 


including lakes from several regions. 
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TABLE 1. Observed average summer chlorophyll] a concentration 
([Chlalsu) and predicted [Chla]lsu based on Dillon and Rigler’s 
(1974) model using spring total phosphorus concentration ([TP]) 
(D&Rsp), Dillon and Rigler’s model using average summer 

[TP] (D&Rsu) and Smith’s (1982) model using average 

summer [TP] and total nitrogen concentration (Cpreéed is the 
predicted [Chla]su and Cobs is the observed [Chlal]lsu). All 


data are from Prepas and Trew (1983). 
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Table 1. Total phosphorus concentration ([TP]), total 
dissolved phosphorus concentration ([TDP]), dissolved oxygen 
concentration ([DO]), and temperature profiles for Nakamun 
Lake (ST = Station). 


DEPTH [TP] [TDP] [Do] TEMPERATURE 
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(Table tt comedy) 


DEPTH [TP] [TDP ] [DO] TEMPERATURE 


(mg/m? ) (mg/m? ) (mg/L) ¢°C) 


0 3.2 62.9 16.4 
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(Table 1 cont'd) 


DEPTH [TP] [TDP ] [Do] TEMPERATURE 
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(Table Tl cont: d) 


DEPTH [TP] [TDP] [DO] TEMPERATURE 
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Table 2. Total phosphorus concentration ([TP]), dissolved 
oxygen concentration ([DO]), and temperature profiles for 
Halfmoon Lake. 


DEPTH [TP] [DO] TEMPERATURE 
(mg/m? ) (mg/L) QC) 
25 May 1982 
0 166.4 Tih, 5 
1 166.4 13.4 11.5 
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7 663.24 0:..0 10.,.6 
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3 55.4 §.a 16.0 
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DEPTH [TP] [DO] TEMPERATURE 
(mg/m ) (mg/L) (2) 
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Table 3. Total phosphorus concentration ([TP]) 
and chlorophyll a concentration ([Chla]) of the 
trophogenic zone, and Secchi disk depth (S.D.) 
in Nakamun Lake. Trophogenic zone defined as 
Z2ucimeseSsbD: 


fChia [TP] Sirs 
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Table 4. Total phosphorus concentration ([TP]) 
and chlorophyll a concentration ([Chla]) of the 
trophogenic zone, and Secchi disk depth (S.D.) 
in Halfmoon Lake. Trophogenic zone defined as 

2, bi meseon. 


[Chla] [TP] Sep 
Date (mg/m? ) (mg/m? ) (m) 
May 27 96.0 PAOOZ O49 
June 41 Zoek 86.9 il 
25 40.8 oe 
Jul 09 54.6 T7837 20 
23 49.9 hea2 .4 
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20 129.4 Uso OFS 
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24 14.0 T3022 23 
Oct 15 Zilee> 166.9 200) 
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Table 5. 


Strata volumes for Nakamun Lake. 


Stratum 


(m) 
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Table 6. Strata volumes for Halfmoon Lake. 


Stratum Volume 


(m) (m? ) 


O= Oro 206173 
1 Bas 850583 
2 aao 324923 
bens} Fis, 297849 
4-4,5 254377 
sees) he 201667 
6- Glico 150087 
Efe h se) 104229 
Om Otel 71874 


Table 7. Total phosphorus runoff loading, aeolean 
loading, loss through the outflow, and net external 
load to Nakamun Lake. 


Runoff Aeolean Outflow NET 

DATE OL -IP TP TP Te 

(kg) (kg) (kg) (kg) 
May 13517 Pa SOE 0.76* uoriz6 eOer a2 
b= 20 Zee 2 0.57% 10.46 SHAS) 
20=20 4.68 0.95% nie OZ Teen! 
EASES h | $307 1.14% 2,49 se 2 
May 3i-Jun 3 0.94 0.69% OFC Lacs 
OT sacs, 0203 0.92% 0.00 Ono 
aat0 0201 0.69 0.00 0.69 
10-14 0.00 0.92% 0.00 0.92 
14-17 0.00 0.69% 0.00 OR69 
UGA 0.00 0.92 0.00 Ue 
21-24 0.00 0.69% O20U 0269 
24-28 0.00 0.92% 0.00 O92 
wun Z28-Jul 1 0.00 0.69% 0.00 Oe) 
Tal) 5 0928 1.24% 0.00 1 5Z 
Sale: O27 0.93% 0.00 eae) 
ape | 2 Org 1.24% 0.00 1.43 
haa 15 a5 25 essa 0.00 6.88 
Posing ns? G..05 0.00 1aG2 
po 1286 O239 0.00 Geo 
22=26 5.54 Oy 512 0.00 6.06 
26729 Zo Ose 0.00 Deny, 
Jthl29-Aug 23 ORs ihe, 1.04 0.00 reo 
Aug 3-6 0.00 0.78 0.00 O76 
ciate, 0.00 2.94 000 2.94 
Set 4 0.00 2.94 6.00 emt fe 
i2aG C200 1.44% 0.00 1.44 
homeo 0.00 1.08% 0.00 Leos 
loa 0.00 1.44% 0.00 1.44 
23526 0.00 he sy, 0.00 toe 
ZO =50 0.00 Zea O00 Zac 
Aug 30-Sep 3 0.00 0.40 O00 0.40 
Sep 3-/ 0.00 io Os00 One 2 
7-14 0.00 Caters 0.00 Orae 
14-21 0.00 O26 0.00 O26 
SBA SAS 0.00 AIS) tie) 0.00 leah! 
Sep 26-O0ct. 5 O00 OneZo 0.00 OAT 
Oct ore 0.00 Oa 0.00 Ori 
2x9 0.00 1.30% 000 Les0 
a2 0.00 eu & 0.00 epee y 
Got /=Now ec 0.00 1.30* 0.00 130 
* Esitmates from average values from Alberta Environment. 
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Table 8. Algal counts from Nakamun Lake, summer of 1982. 


Cyanophytes Chlorophytes Cryptophytes Diatoms 


Gells/m. Cells/mL Cells/mL Cells/mL 

Date ya 0- EPO ke Xe ropes 
May 17 0.00 2.60 0.00 765.00 
25) O00 Sion) aden) JOS. 00 
31 2.80 1.90 Se) CLR SHG) 
June 27 37.40 5, 60 10 Soo 
14 2297700 9.70 BT T6) 0.00 
2a B25n06 Ha 200 ohn JAD) O00 
28 305.00 9.40 O00 0.00 
te 222700 cpa ic 8 S90 HeoU 
19 250700 3. 90 1.50 Oe10 
PAS) 5747..00 3.90 e550 0.00 
Aug 2 PSA O wa ORO 0.00 0.00 0.00 
9 S200 2.40 0.00 OF 0 
16 340.00 0.80 O00 0.00 
Ze 242.00 550 G00 0.00 
30 PRI EERONe) O00 0.00 0.00 
Sep ZOO U0 0.00 0200 2.30 
14 198.00 0.00 OE, 0.0 0.00 
20 VOR S00 Zoe 9250 G00 
OGraao WiPet, Oe 0.00 ES290 0.00 
12 152.90 0.00 7.90 0.00 
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